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Abstract: It is shown that lattice kinetic theory based on short-lived quasiparticles proves 
very effective in simulating the complex dynamics of strongly interacting fluids (SIP). In 
particular, it is pointed out that the shear viscosity of lattice fluids is the sum of two 
contributions, one due to the usual interactions between particles (collision viscosity) and 
the other due to the interaction with the discrete lattice (propagation viscosity). Since the 
latter is negative, the sum may turn out to be orders of magnitude smaller than each of the 
two contributions separately, thus providing a mechanism to access SIP regimes at ordinary 
values of the collisional viscosity. This concept, as applied to quantum superfluids in 
one-dimensional optical lattices, is shown to reproduce shear viscosities consistent with the 
AdS-CPT holographic bound on the viscosity/entropy ratio. This shows that lattice kinetic 
theory continues to hold for strongly coupled hydrodynamic regimes where continuum 
kinetic theory may no longer be applicable. 
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1. Introduction 

The study of transport properties of strongly interacting fluids (SIP) has gained central stage in modern 
condensed matter research, with many fascinating connections with quantum-relativistic hydrodynamics, 
high-energy physics and string theory [1]. By definition, SIP are moving states of matter in which the 
interactions are strong to the point of preventing the microscopic degrees of freedom from propagating 
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freely over signifieant distanees as compared to the range of their interactions. More precisely, the 
mean-free path becomes smaller than the interaction range, and, eventually, the De Broglie length. 
Remarkable examples in point are found mostly in the framework of quantum fluids, namely quark-gluon 
plasmas [2], electrons in graphene [3] and Bose-Einstein condensates in optical lattices [4]. Note that 
these extreme and exotic states of matter cover a breathtaking range of densities and temperatures, 
from one particle/fm^ and trillions Kelvin degrees in ultra hot-dense quark-gluon plasmas, to about 
10^^ particles/cm^ and tens of nanokelvins in ultra-cold Bose-Einstein condensates. Thus the SIF 
regimes show an impressive degree of universality, whose understanding stands as a great challenge 
at the crossroad of statistical fluid mechanics, condensed matter and high-energy physics. Under SIF 
conditions, most of the familiar and powerful notions of kinetic theory call for a profound revision. In 
particular, it is argued that not only Boltzmann’s kinetic theory, but the very notion of quasi-particles as 
weakly-interacting collective degrees of freedom, would loose meaning due their inconspicuously short 
lifetime [5]. Hence, new ideas and methods, both analytical and numerical are in great demand. 

Among the analytical ones, a most prominent role has been taken by the AdS-CET (Anti de Sitter 
- Conformal Eield Theory) duality between gravity in (d -I- 1) -dimensions and conformal field theory 
in d-dimensions [6] . The beauty, and practical import, of the dual-holographic picture is that one can 
solve an otherwise intractable strongly coupled CET by dealing with its weakly coupled gravitational 
analogue, an approach sometimes known as holographic principle. One of the major outcomes of the 
holographic approach is the so called minimum-viscosity bound [7] 


rj ^ 1 h 
s ~ An ks 


( 1 ) 


where rj is the shear viscosity and s the entropy per unit volume. It has been found that many 
quantum relativistic fluids, such as quark-gluon plasmas and electrons in graphene, come much closer to 
matching the above bound than any previously known fluid, including superfluid Helium-III. Moreover, 
the holographic AdS-CET picture shows that hydrodynamics continues to apply even though ki netic 
theory does (may) not. Besides the conceptual challenge of developing statistical mechanics without 
quasiparticles, the ultra-low values of p/s open up the exciting possibility of observing dynamical 
instabilities and ensuing (pre) quantum turbulence phenomena in these extreme states of matter [8]. 

It is hereby maintained that a proper lattice formulation of effective kinetic theory can effectively 
describe the dynamics of SIE’s by making use of very short-lived quasiparticles. The two main 
conceptual ingredients are as follows: 1) Top-down approach: design kinetic equation solely based 
on symmetry and conservation properties of the macroscopic field-theory (hydrodynamics), including 
entropic constraints. This just the reverse of the canonical route of deriving kinetic theory from 
first-principles, i.e. by coarse-graining underlying microscopic dynamics (which is the process 
generating quasi-particles). 2) The existence of negative propagation viscosity, as an emergent property 
due to the interaction of the free-streaming quasiparticles with the discrete lattice. Item 2) is crucial to 
sustain ephemeral quasiparticles, with lifetimes as small as one millionth of the free-flight time, hence to 
achieve ultra-low viscosities, well below the natural kinematic lattice viscosity ul = Ax^/Af, Ax and 
At being the lattice spacing and hopping time, respectively. 
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2. Results and Discussion 

The typical lattice kinetic equation reads as follows [9] : 

-/i(f - - At) =-^[/i -- At), (2) 

where /i(x, t) is the probability of finding a particle at position x in the lattice at time t with discrete 
velocity v = Ci. The left hand side represents the free-streaming from site to site while the rhs encodes 
collisional relaxation to a local equilibrium distribution on a time scale r (see Fig. 1). The discrete 
Boltzmann distributions (“populations”) are quintessential quasiparticles, as they encode the dynamics 
of a collection of molecules in a lattice of size Ax^ (in three dimensions), moving along the direction 
Ci- The local equilibrium is a universal function (e.g. Maxwell-Boltzmann, Bose-Einstein, Fermi-Dirac, 
etc) of the locally conserved quantities (hydrodynamic modes), such as flow density p = m fi 
current J = pu = m Aft, and dictates the structure of the inviscid macroscopic equations (Euler). 
Here m denotes a characteristic mass. The relaxation time, on the other hand, corresponds to the rate at 
which this equilibrium is reached, and dictates the fluid transport coefficients. 

In the Boltzmann kinetic theory, in the strong-coupling regime r —)■ 0, viscosity is directly 
proportional to this relaxation time, v ~ c^r, Cg = y/kBT/m being the thermal speed [10]. In 
the lattice, however, this expression receives a crucial contribution from the discrete free-streaming. 
Indeed, by performing the appropriate Chapman-Enskog asymptotic expansion in the Knudsen number 
Kn = Cst/L, L being a typical hydro-scale, one obtains: 

^ = cliT-At/2) , (3) 

where the negative term at the rhs is the contribution of the free-streaming. Given the minus sign, it is 
immediately seen that ultra-low viscosities of order can be achieved by near-matching the free-flight 
and collisional time, i.e. by choosing r = Af(l -f e)/2. In lattice units, one can write Eq. (3) as 

o = ClMl/uj-l/2) , (4) 

where ul = is the natural lattice viscosity and we have set cu = Af/r. In the above Cg is the sound 
speed in lattice units, typically l/\/3 in most common lattices. Crucial for the correct recovery of the 
hydrodynamic equations is that the lattice exhibits enough symmetry to ensure the proper conservation 
laws as well as rotational invariance [11]. 

To further appreciate the point of ultra-low viscosity, let us now recast the kinetic equation (2) in the 
most compact stream-collide form: 

fi{x,t) = f-{x - CiAt,t - At) , (5) 

where 

/• = (!- uj)fi + , ( 6 ) 

is the post-collisional distribution (see Fig. 1). The latter expression highlights three distinguished 
updates: i) u = 0, fi = fi. The particles are fully uncoupled (no collision) perform ballistic motion 
and never equilibrate, corresponding to infinite diffusivity. ii) cu = 1, // = ff^'. The particles reach 
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Figure 1. Schematics of the LB scheme: the quasi-partiele distributions hop along the 
links defined by the discrete velocities and, onee on the same site, they scatter into the 
post-collisional state. One- (top) and two- (bottom) dimensional lattices. Here Ax is the 
size of the lattiee. 
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equilibrium in a single time-step. According to Eq. (4), this gives a viscosity v = z/iC^/2, i.e. 
comparable with the lattice viscosity ol. iii) u = 2, /' = fi — 2/”®'^; Collisions send the populations 
exactly “on the other side”, to the mirror state f* = defined by complete reversal of the 

non-equilibrium component = /j — Based on Eq. (4), this yields formally zero viscosity, i.e. 
infinitely strong coupling. 

Hence, by choosing r = Af/2, the quasi-particles have literally zero life-time, even though both 
free-streaming and collisional time-scales are finite. Ear from being a mere mathematical nicety, this 
regime is crucial to the operation of the EB scheme in the very-low viscous regime relevant, say, to fluid 
turbulence. Indeed, given the current computer resolutions, iii) is precisely the regime that needs to be 
approached in order to simulate strongly coupled classical fluids, such as turbulent flows at Reynolds 
number above a few thousands. To clarify the point, let us recall that the strength of fluid turbulence is 
measured by the Reynolds number. Re = UL/o, where U and L are macroscopic velocity and length 
scales, respectively. In lattice Boltzmann units. Ax = At = 1, this gives Re = MaN/un,, where 
Ma = U/cs is the thermal Mach Number, N = L/Ax is the number of lattice sites per linear dimension 
and oi}, = v/vl = (l/cu — 1/2), is the viscosity in lattice Boltzmann units. Current computers allow at 
best N ~ 10"^, so that with Mach numbers of order 1, reaching i?e = 10^ (air flowing around a standard 
car), requires uib ~ 10“^, i.e. three orders of magnitude below the natural lattice viscosity i>l- This 
simple example highlights the necessity of operating with very short lived quasiparticles, where short 
means very small as compared to the physical free-streaming and collisional time scales. At and r, both 
0(1) in lattice Boltzmann units. 

These quasi-particles function pretty well in the discrete world they live in, if only very shortly. In fact, 
by promoting the relaxation parameter oj to the status of a dynamic field responding self-consistently to 
the local constraints of the second principle (Boltzmann’s H-theorem), the lifetime of these “ephemeral” 
quasi-particles can be made as small as one millionth of their natural value, i.e u; ~ 2(1-10-®) [12-14]. 
In order to appreciate why this is remarkable, let us consider the lattice kinetic update in compact form, 
Eq. (5). The realisability constraint imposes that the discrete populations be non-negative, hence, starting 
with a non-negative pair (/j, //'^), the collisional update should return a non-negative //. A moment 
thought shows that this is indeed the case in the safe under-relaxation regime 0 < ce < 1. However, 
we have just shown that such regime does not provide access to the large Reynolds numbers typical of 
most turbulent flows. Hence, the EB update needs to operate deeply into the potentially unprotected 
over-relaxation regime, 1 < ce < 2. The remarkable fact is that this is indeed possible, thanks to the 
existence of a lattice analogue of the H-theorem [12-14]. Note that the regime 1 < cj < 2 is called 
over-relaxation because the time step is larger than the relaxation time, r < At, and consequently, 
the post-collisional distribution is no longer bounded between the pre-collisional one and the local 
equilibrium, which is a potential threaten to positive-definiteness. 

2.1. Propagation viscosity in Luttinger liquids 

Eet us apply this concept to Euttinger liquids. By dealing with a unitary Eermi gas, one should 
distinguish two different regimes, weak and strong coupling. Since the present work is focused on 
strongly interacting fluids, we consider the regime of large values of the Eieb-Einiger parameter, 7 = 


k = {B! jt)l 

<—> 



Ag = 4/ 


Figure 2. Representation of the Luttinger liquid. Shown are the mean-free-path A, the De 
Broglie length As, the lattiee size Ax = I, and the effeetive mean-free-path predieted by 
the lattiee Ag//. Note that, due to the negative contribution {—1/2), lattice kinetic fluids can 
attain ultra-low viscosity even though the standard collisonal viscosity is of the same order 
of the lattice viscosity vl- 


mg/{h^n) 3> 1, where m is the atomic mass, n the atomic density, and g the coupling constant [15]. 
It proves expedient to define the dimensionless quantity K = hmi/(mCs). Note that 7 and K relate to 

1 /*A 

each other through an expression that, for 1 < 7 < 10, takes the form K ~ vr [7 — (l/27r)7^/^] 

[15]. 

In the strongly interacting regime and for weak external potentials, a quantum phase transition from 
Mott-insulator to superfluid occurs for 7 > 7 c and = 2 [15,16]. 

Let us next define, from Eq. (3), z/g = c^r and Vp = At/2, denoting “collision” and “propagation” 

viscosities, respectively. In order to estimate the propagation viscosity for a SIF at quantum criticality, 
we recast Eq. (3) in the form: 

i/p/i/c = At/ 2 r . (7) 

Assuming that At = l/c, where I is the lattice periodicity, with c the light speed, and r = A/c^, with A 
the characteristic mean-free path, the above relation reads: 

z^p/z/c = (l/2)(cg/c)(//A) . ( 8 ) 


In the sequel, we shall take c ~ Cg. Since at criticality, Kc = 2, taking into account that the density 
n ~ 1/1 (commensurate density [15]), it follows that Xb = 'll, where Xb is the De Broglie length. 
Inserting this into the viscosity ratio, we obtain 


z/„ 


crit 


1 Xb 

sT' 


(9) 


Note that SIFs regimes are characterized by the condition A < A^. To determine the ratio A^/A of the 
lattice fluid, we follow the same procedure adopted in Ref. [18] to obtain the AdS-CFT viscosity bound. 
Based on the Heisenberg principle, /c^Tr > h/2, so that r = A/c = fiB/{2kBT), B > 1 being a 


constant. 
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Consequently, 

A = ^\b . (10) 

47r 

The relations between the lattiee size I, the mean-free path A, the De Broglie length As and the predieted 
lattice mean-free path Ae//, can be appreciated from Fig. 2. 

Note that the constant B is related to the AdS-CFT bound by [4], 


r] B h 

s dvr kB 


( 11 ) 


Inserting Eq. (10) into Eq. (9), we obtain 


Urt 


crit 


TT 

w 


( 12 ) 


This result calls for a number of comments. First, we note that the condition of non-negative viscosity, 
t'p/t'c < 1, implies B > n 12 ~ 1.6, which corresponds to a bound rf/s ~ O.lS^/fcs. In this case, 
the lattice quasi-particles can model even more "perfect" fluids than the Luttinger liquid. On the other 
hand, if we assume that r = At, which means that the quasi-particles thermalise in a single time step, 
we obtain B = tt, yielding a bound rj/s ~ 0.25 h/kB, which is very close to the value measured 
in experiments for two and three dimensional Fermi gases at unitarity [17]. Note that in order to get 
Eq. (12) we have not only used lattice kinetic theory but also the Heisenberg principle, which shows that 
we still need quantum mechanics contributions for this particular example. However, this does not imply 
that lattice kinetic theory cannot handle quantum systems by itself. In fact, a quantum lattice Boltzmann 
model exists since long [20] to solve the Dirac equation. It is based on the same stream-collide paradigm 
of the classical lattice Boltzmann equation, although with a skew-symmetric collision matrix, where r is 
related to h via r = h/mc^. 


3. Conclusions 

In summary, one can conclude that lattice kinetic schemes endowed with an H-theorem support 
ultrashort-lived computational quasiparticles, down to viscosities several orders of magnitude below 
their natural scale (the lattice viscosity). This is a strict consequence of the fact that, due the discrete 
nature of the lattice, free-streaming also contributes a viscosity and, most importantly, a negative one. 
It is as if, in the process of hopping from site to site, the lattice quasiparticles would “surrender” part 
of their collisional viscosity to the lattice itself, and precisely in proportion of half of the flight-time. 
This is plausible, since negative viscosity corresponds to dynamic instability, namely the system moving 
away from local equilibrium, a sort of backward move in time (rejuvenation), as opposed to the standard 
forward move due to collisions. It is also plausible to speculate that this ultra-low viscosity could be 
regarded as an effect of lattice criticality, in the sense that the two competing processes, collisions and 
streaming, come to a near exact balance through the interaction of the quasiparticles with the lattice. 

This process follows in the footsteps of similar phenomena in condensed matter: electrons in graphene 
behave like effective near mass-free excitations precisely because, due to the special symmetry of the 
honeycomb lattice, the electrons near the Dirac point “release” most of their mass to the lattice [19]. 
Since mass dictates the (Compton) collision frequency, uOc = me?/%, and the collisional viscosity is 
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inversely proportional to ujc, it is elear that negative viseosity and near mass-free eleetrons in graphene 
belong somehow to the same family of lattiee-indueed phenomena. In the eontinuum, At/r —?■ 0, the 
negative viseosity is totally negligible, but in the diserete world it beeomes erueial to attain ultra-low 
kinematie viseosities. 

By using lattiee kinetie theory, one ean indeed aehieve zero viseosity, and the same ean be done in 
the eontinuum theories, sueh as eontinuum Boltzmann and Navier-Stokes equations. However, when 
the physieal system possesses intrinsie diserete properties, e.g. quantum liquids in optieal lattiees, 
the viseosity ean take very small values due to the presenee of the negative viseosity indueed by the 
lattiee. For these partieular systems, lattiee kinetie theory may deseribe better the transport properties 
than eontinuum kinetie theory, as is the ease of the lower bound in the shear viseosity of quantum liquids, 
whieh is overestimated by the eontinuum Boltzmann equation. On the other hand, in the elassieal 
eontext, a perfeet fluid has zero viseosity, whieh eorresponds in our ease to ce = 2. Thus, in the 
eontinuum limit, the lattiee kinetie theory beeomes just a numerieal tool and reeovers the Euler equations. 
However, in strongly eoupled fluids, whieh are the focus of this paper and usually are governed by 
quantum mechanics, fluids do not have zero but a finite value of viscosity. Here we show that considering 
cj = 1 one can get instantaneous relaxation and at the same time a lower bound for the viscosity of the 
physical system. 

A few concluding remarks regarding strongly coupled fluids are in order. One may argue that most 
of these fluids are quantal, hence not captured by classical lattice kinetic models, for instance, the LB 
schemes discussed in this work. To this regard, it is worth pointing out that the lattice Boltzmann method 
has been already used to study two-dimensional Fermi gas at unitarity [? ], providing a procedure to 
calculate transport coefficients precisely. Furthermore, LB versions for quantum wave functions are 
available since long, and they are based exactly on the same stream-collide paradigm, if only with a 
suitable scattering matrix for collisions [20]. Incidentally, such quantum LB versions provide exact 
realisations of the Dirac equation in one spatial dimension and can be extended to d > 1 by proper 
operator splitting [21]. Finally, we wish to point out that matrix extensions of the quantum LB scheme 
have been recently connected with the Hubbard model for strongly correlated quantum fluids [22], 
thereby lending further weight to the main idea proposed in this paper, namely that lattice kinetic theory 
may provide a new angle of investigation for strongly-interacting fluid regimes where continuum kinetic 
theory may no longer hold. 
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